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The  Collinear  Crack  Problem  in  an 
Inhomogeneous  Orthotropic  Medium 

by 

Murat  Ozturk  and  Fazil  Erdogan 
Department  of  Mechanical  engineering  and  Mechanics 
Lehigh  University,  Bethlehem,  PA  18015 

Abstract 


The  collinear  crack  problem  in  an  inhomogeneous  orthotropic  medium  is 
considered  under  Model  plane  strain  or  plane  stress  loading  conditions.  It  is  shown 
that  by  introducing  certain  averaged  orthotropy  parameters,  aside  from  a  scaling 
parameter  the  results  become  only  weakly  dependent  on  the  orthotropy  constants.  Die 
main  results  of  the  study  consist  of  the  stress  intensity  factors  at  various  crack  tips  as 
influenced  by  tlie  material  inhomogeneity  parameter  and  by  the  relative  size  and 
position  of  the  secondary  cracks  with  respect  to  the  dominant  crack.  Also  considered 
is  the  crack/contact  problem  for  the  graded  medium  subjected  to  remote  tension  and 
bending  through  fixed  grips.  It  is  shown  that  the  crack  surface  contact  on  the 
compression  side  of  the  loading  has  a  magnifying  effect  on  the  stress  intensity  factor 
on  the  tension  side. 


1.  Introduction 

In  recent  years  there  has  been  considerable  interest  in  grading  the  thermomechanical 
properties  of  composites  as  a  new  tool  in  designing  materials  for  specific  applications  ( see 
Yamanouchi  et  al.  1990,  Holt  et  al.  1993  and  Ilschner  and  Cherradi  1995  for  review, 
applications  and  extensive  references).  Up  to  now  most  of  the  work  in  the  field  has  been 
on  metal/ceramic  particulate  composites.  By  proper  selection  of  the  constituent  materials, 
compositional  grading  and  the  processing  technique,  the  concept  may  be  used  to  develop 
new  materials  having  such  highly  desirable  and  seemingly  irreconcilable  properties  as  high 
heat,  corrosion  and  wear  resistance,  high  strength  and  high  toughness  in  the  same  material 
system.  The  technique  may  also  be  used  to  process  interfacial  regions  to  improve  bonding 
(Kurihara  et  al.  1990)  and  to  reduce  the  magnitude  of  residual  and  thermal  stresses  in 
bonded  dissimilar  materials  (Lee  and  Erdogan  1995,  Lee  and  Erdogan  1997). 

In  the  past  the  fracture  mechanics  studies  of  graded  materials  were  concerned 
primarily  with  the  influence  of  material  inhomogeneity  constants,  and  certain 
dimensionless  length  parameters  on  the  fracture  behavior  of  components  containing  a 
single  dominant  crack.  For  example,  the  model  and  the  mixed  mode  crack  problems  for 
an  infinite  isotropic  inhomogeneous  medium  were  studied  by  Delale  and  Erdogan  (1983) 
and  Konda  and  Erdogan  (1994),  respectively.  The  influence  of  the  length  parameters  on 
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the  stress  intensity  factors  and  the  strain  energy  release  rates  in  isotropic  graded  layers 
undergoing  spallation  and  surface  cracking  was  investigated  by  Chen  and  Erdogan  (1994) 
and  Erdogan  and  Wu  (1996),  respectively.  However,  because  of  the  nature  of  the 
techniques  used  in  processing,  the  graded  materials  are  seldom  isotropic.  Of  the  two  most 
commonly  used  processing  techniques,  generally  thermal  spray  would  give  a  lamellar 
(Sampath  et  al.  1995)  and  the  electron  beam  physical  vapor  deposition  a  columnar 
structure  (Kaysser  and  Ilschner  1995).  An  appropriate  model  for  such  graded  materials 
would  be  an  orthotropic  inhomogeneous  continuum  (Ozturk  and  Erdogan  1997).  One 
additional  factor  that  needs  to  be  considered  in  studying  the  basic  crack  problem  in  graded 
materials  is  the  so-called  multiple-site  cracking  or  the  crack  interaction.  In  homogeneous 
materials,  as  intuitively  expected,  the  crack  interaction  would  invariably  increase  the 
relative  magnitude  of  stress  intensity  factors.  In  graded  materials,  however,  such  a 
categorical  statement  is,  generally  not  possible  and  the  result  may  depend  on  the  relative 
positions  of  the  multiple  cracks  with  respect  to  the  direction  of  the  material  property 
variation  as  well  as  on  the  relative  dimensions. 

The  main  objective  of  this  article  is  to  study  the  influence  of  the  relative  size  and 
location  of  a  secondary  crack  on  the  stress  intensity  factors  in  an  orthotropic 
inhomogeneous  medium  containing  a  dominant  crack.  It  is  assumed  that  the  cracks  arc 
collinear  and  are  located  in  a  plane  parallel  to  the  direction  of  material  property  variation. 
It  is  also  assumed  that  the  problem  is  one  of  plane  strain  or  generalized  plane  stress  and 
the  solution  of  the  elasticity  problem  in  the  absence  of  crack  is  known.  Thus,  the  plane  of 
the  cracks  is  one  of  symmetry  and  the  problem  is  a  model  crack  problem  in  which  the 
known  crack  surface  tractions  are  the  only  external  loads.  The  problem  is  formulated  for 
arbitrary  crack  surface  tractions  and  the  results  are  given  for  some  simple  loading 
conditions.  The  stress  intensity  factors  for  more  complex  loadings  may  then  be  obtained  by 
superposition.  Clearly,  the  result  would  be  valid  only  if  the  model  stress  intensity  factors 
kx  at  all  crack  tips  are  positive.  A  negative  ki  implies  crack  closure  in  which  case  the 
problem  becomes  nonlinear  with  the  size  of  the  contact  zone  being  an  additional  unknown. 
This  problem  is  solved  by  assuming  that  the  graded  medium  contains  a  single  crack  and  is 
subjected  to  tension  and  bending  away  from  the  crack  region  through  rigid  grips  (  see 
Appendix  A  for  the  solution  in  a  homogeneous  isotropic  medium  which  is  obtained  in 
closed  form). 

2.  Formulation  of  the  Elasticity  Problem 

Let  Exx,  E22,  Gx2  and  Vij,  (i,j  =  1, 2, 3)  be  the  engineering  elastic  parameters  for  an 
orthotropic  inhomogeneous  plane  xi,  X2.  To  replace  them  we  introduce  the  following  four 
parameters  (Krenk  1979,  Cinar  and  Erdogan  1983): 

E  V  E 

E  =  y/ExxE22,  V  =  ^  ^  ^  2^  ~  ^ 


for  plane  stress  and 
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_ EIIE22 _ 

(1  —  I^i3^'3l)(l  ~  ^23^32) 


1/2 


{Ui2  +  ^13  ^32)  (^21  +  ^23  ^31 ) 

(1  —  Ui3  1^31  )(1  —  1^23  1^32) 


1/2 


,l^  = 


? 


^11  (1  ~  1^231^32) 
E22  (1  ~  f^l3Z^3l)’ 


K  = 


E 

2G\2 


(2) 


for  plane  strain  conditions.  In  graded  materials  generally  E,  u,  k  and  So  would  be 
functions  of  xi,X2.  However,  previously  it  was  shown  that  the  results  are  only  weakly 
dependent  on  ly  (Delale  and  Erdogan  1983,  Ozturk  and  Erdogan  1997)  and  it  may  be 
assumed  that  En,  E22  and  Gu  vary  proportionately.  Consequently,  the  parameters  k  and 
So  may  be  assumed  to  be  constant  and  the  material  inhomogeneity  may  be  represented  by 
the  function  E{xi,X2)  only.  Note  that  5o  =  1  and  k  =  1  correspond  to  the  isotropic 
medium  and  for  k  <  -1  the  problem  has  no  feasible  solution  (Erdogan  and  Wu  1993).  By 
introducing  So  as  a  scaling  parameter,  in  the  usual  notation  we  define 

X  =  xi/ y  =  ^2, 

u{x,y)  =  Ul{xl,X2)^/^,  v{x,y)  =  U2{x\,X2)/ 

(^xx{x,y)  =  cru{^i,X2)/So,  (^yy{x,y)  =  So  a22(xi ,  X2) , 

(^xy{x,y)  =  ai2{Xl,X2).  (3) 

Furthermore,  by  assuming  that  the  stiffhess  E  varies  in  xi  direction  only  and  in  the  crack 
region  may  be  approximated  by  (Fig.  1) 

EixuX2)  =  Eixi)  =  EQe°‘^^  =  E{x)  =  Eoe^\  7  =  a\/^,  (4) 


the  equilibrium  equations  may  be  expressed  as 


d'^u  d^u  d'^v 
dy"^  ^  dx^  ^  ^  dxdy 


+  At 


du  dv\ 
dx  ^  dy  ) 


=  0, 


d^v  -  d^v  „  d'^u  f  du  dv\ 
dx^  '  dy"^  ^ dxdy  \ dx) 


=  0, 


(5) 


where 


1-1^2  > 

Due  to  symmetry  by  considering  y  >  0  half  plane  only,  from  (5)  it  may  be  shown  that 


(6) 


1  r®®  ^ 

»)  =  -y  ^  ^  dk, 

-I  poo  4 

V)  =  ^J  E  0  <  »  <  0°. 


(7) 
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(^xy{x,y) 


(8) 


where 


Ai  =  —  A3  =  I  -  (C^  +  2«7?)  +  -  \/ (C^  +  2«7;)^  —  4.7]'^  I  , 

A2  =  -  A4  =  ||(C^  +  2Ky)  -  (C^  +  2KT]f  -  477^1  ,  (9) 

=  i>'7^,  ri  =  k‘^  +  i'yk,  (10) 


and 


(e/32 A;  -  . 


J  =  3,4. 


(11) 


The  unknown  functions  and  are  determined  from  the  following  boundary 
conditions 


(7’i2(a:i,0)  =  0,  -00  <  a;i  <  00.  (12) 

n 

(^22(^1, +0)  =  p{xi),  xi  e  L,  L  =  Y^Lk,  Lk  =  {akbk), 

1 

^2(3^1? +0)  =  y^t'(a;, +0)  =  0,  xiEL',  L  +  L' =  {—00,00).  (13) 

Here,  it  is  assumed  that  the  crack  surface  traction  p{xi)  is  a  known  function  and  the 
medium  contains  n  collinear  cracks  along  X2  =  0  (Fig.  1). 

By  defining  the  following  new  unknown  function 

- — «2(a:i,  +0)  =  gixi),  -00  <  Xi  <  00,  (14) 

UXi 

it  can  be  shown  that  the  boundary  conditions  reduce  to  (Ozturk  and  Erdogan  1997) 


Xi 


+  7^2 (ii  ~  ^i) 


g{ti)dti 


(si  +  S2)5o  _axi 


Eq 


e  "^'p(a;i),  xi  e  L, 


(15) 
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g{xi)dxi 


=  0, 


(16) 


where 


Si  =  a/aC  +  ACi,  S2  =  y'/C  -  Kl,  -  1, 

■^2(^1  ■”  2^1)  =  ((^1  ~  2;i)/\/^) 


(17) 


i?i(s-a;)=  /  2  [3?(/i(/:))cos(A:(s  -  a:))  -  Q‘(/i(A:))sin(/i:(s  -  a;))  dk,  (18) 


^(A:)  = 


— e 


(si  +S2)(fc  +  i7)  _  ' 
Ai  +  As 

2{k  +  i^) 


-  1 


/c  7^  1 


K  =  1. 


(19) 


In  real  materials  -1  <  /«  <  ooand  Si  +  ss  is  always  real.  Also,  for  7  =  0,  Ai  =  ksi, 
A2  =  ks2,  h{k)  =  0,  i?2  =  0  and  (15)  would  reduce  to  the  integral  equation  for  a 
homogeneous  orthotropic  medium. 

It  is  interesting  to  note  that  in  the  special  case  of  isotropic  inhomogeneous  materials 
the  integral  in  (18)  can  be  evaluated  in  closed  form  giving  (Appendix  B) 

R2{h  -  Xy)  =  +  Ko{uo\^\)^  - 


^=^(ti-xr). 


(20) 


where  Kq  and  Ki  are  modified  Bessel  functions.  Also  note  that  the  kernel  R2  depends  on 
the  inhomogeneity  constant  a  and  the  elastic  constants  u  and  k  but  not  on  Eq  and  60  • 
Similarly  si  and  ss  depend  on  the  shear  parameter  k  only.  It  may  then  be  concluded  that  in 
(15)  Eo  and  60  are  simply  scaling  constants  and  would  have  no  influence  on  the  stress 
intensity  factors  which  may  be  defined  by  and,  after  solving  (15),  evaluated  from 


ki{bk)=  lim  \/2(x7^^^^a22(a:i,+0) 

^  h  +0 

=  -lim  -  xi)  g{xi), 

ii-6*-o(si +52)^0 


ki  {flk )  =  lim  y/2{ak  -  Xi)  <752 (a:i ,  +0) 

Xi  — >  ak  —0 


\/2(2;i  -  afc)  ^(xi),  fc  =  1, 

xi  — ^  ak  +0  \^Si  +  52  jOq 


(21) 
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(22) 


Also,  from  (13)  and  (14)  the  crack  surface  displacements  may  be  obtained  as 

PXl 

^2(^1 )  ~b^)  ~  )*^^i »  ^  ^1  ^  3  1)  •  •  • ) 


3.  On  The  Solution  Of  The  Integral  Equation 

Since  the  closed  form  solution  is  not  available  the  integral  equation  (15)  must  be 
solved  numerically.  A  simple  and  very  effective  way  to  do  this  would  be  to  reduce  (15)  to 
a  system  of  n  integral  equations  each  having  the  support  (—1,  1)  by  defining 


2^1  —  {bk  +  ak) 

^  {bk  —  o^k)  ’ 

^1  ^  {^k> 

-1  <q<l. 

23^1  (fe/c  H"  Cl/j;) 

ak  <  xi  <  bk, 

—  1  <  r  <  1, 

1 

1 

<t>k{q)  =  ak<ti<hk,  -1  <  ?  <  1, 

fk{r)  =  p{xi)/Pq,  ak  <xi  <  bk,  -1  <  r  <  1, 

RikMq  -  Zik))  =  R2{h  -  xi)/ck,  ai<xi<bi,  ak<ti<bk, 

-l<q<l, 


Cfc  =  2^bk  —  ak),  dk  =  2 

Zik  ^ —r  + —{di  -  dk),  {i,k)  =  -1  <  r  <  1,  (23) 

Ck  Ck 


where  Po  which  has  the  dimension  of  stress  is  a  normalization  constant.  The  equations 
(15)  and  (16)  may  then  be  expressed  as 


1 

TT 


“f*  OiRii 

q  —  r 


Mq)  dq 


+ 


1 

.q  -  Zik 


“b  ^kR-ik 


(l>k{q)dq  = 


(si  +  52)^oPo  r  /  n 
£;oexp(a;(cir  +  di)) 


-1  <  r  <  1,  i  =  (24) 


/: 


(t>k{q)dq  =  Q,  k  =  l,---,n. 


(25) 


The  system  of  integral  equations  may  be  solved  in  a  standard  manner  by  expressing 
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0.(9) 


(si  +  S2)6oPo  1 

£;oexp(Q;4) 


■^kmTm(,Q\ 


—  l<q<l,  k  = 


(26) 


where  Tm{q)  is  the  Chebyshev  polynomial  of  the  first  kind  and  Akm  are  unknown 
dimensionless  constants.  From  (25)  and  othogonality  conditions  of  T^'s  it  follows  that 
AkQ  =  0.  The  coefficients  Akm  are  determined  by  substituting  from  (26)  into  (24)  and  by 
using  the  method  of  reduction  (Kantorovich  and  Krylov  1958,  Erdogan  1978).  The  stress 
intensity  factors  may  then  be  obtained  from  (21)  and  (26)  as  follows; 

00 

h{hk)  =  Akm, 

Tn~l 


h{ak)  =  PoV^e-“‘=*J](-l)™Afe^,  k  =  l,--;n.  (27a, b) 

m=l 


4.  Results  and  Discussion 

The  main  interest  in  this  study  is  in  evaluating  the  influence  of  relative  position,  size 
and  distance  of  a  "small"  crack  lying  in  the  plane  of  a  dominant  crack  on  the  stress 
intensity  factors  in  an  inhomogeneous  orthotropic  medium.  With  the  additional  knowledge 
of  subcritical  crack  growth  characteristics  of  the  material,  one  may  then  be  able  to 
determine,  for  example,  fatigue  or  corrosion  crack  growth  rate  at  each  crack  tip.  The 
problem  is  solved  under  two  sets  of  loading  conditions.  In  the  first,  it  is  assumed  that  the 
crack  surface  tractions  are  given  by  (see  13) 

p(i,)=  (28) 

where  po,Pi,P2  and  ps  are  the  measure  of  the  magnitude  of  external  loads  and  are  known 
constants  and  a  is  a  normalizing  length  parameter  (usually  the  half  length  of  the  dominant 
crack  ).  Since  the  problem  is  linear,  these  results  may  be  helpful  to  find  an  approximate 
solution  to  the  crack  problem  for  a  given  specific  loading.  The  second  external  load 
considered  is  a  remote  displacement  loading  of  the  form 

£22(3)1, 0)2)  =  Co  +  ei(a;i/a)  (29) 

giving 

p(x, )  =  -  (eo  +  €1  (^) ) ,  3:,  e  L.  (30) 

Summary  of  the  results  for  an  isotropic  inhomogeneous  medium  containing  a  single 
crack  along  X2  ^  0,  -a  <  xi  <  a  is  shown  in  Tables  1-3.  These  results  are  both  accurate 
and  comprehensive.  Tables  4-6  show  the  comparison  of  the  stress  intensity  factors 
obtained  in  this  study  with  that  found  by  Delale  and  Erdogan  (1983).  From  the  solution 
given  in  the  previous  section  the  crack  surface  displacement  may  easily  be  expressed  as 
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U2{xi,  +0) 


(31) 


where  112(2:1,  +0)  =  v{x,  0), 

Vo  =  2aPo{l-u'^)/Eo  (32) 

for  isotropic  and 

Vo  =  (si  S2)SoPoCi/ Eq  (33) 

for  orthotropic  inhomogeneous  materials  (Eq  =  Eq  for  plane  stress  and 
Eo  =  Eo/{l  —  for  plane  strain  ).  Some  sample  results  for  the  crack  surface 
displacement  are  shown  in  Figure  2. 

The  calculated  plane  strain  results  for  the  normalized  stress  intensity  factors  at  various 
crack  tips  in  an  isotropic  inhomogeneous  medium  containing  two  unequal  cracks  under 
fixed  grip  (£22  =  ^o)  or  fixed  load  (cr j,j, (a: i,0)  =  -po)  conditions  are  given  in  Figures  3- 
16.  Some  of  the  limiting  cases  of  these  results  for  (ci  +  C2)/d  0  (the  uncoupled  case) 

and  (ci  +  C2)/d  — » 1  ( the  case  of  a  single  crack  of  length  2(ci  +  C2))  may  be  found  in 
Table  7.  One  may  note  that  as  (ci  +  C2)/d  — >  1  or  a2  — >  61,  ki{bi)  00,  ki{a2)  00 
and  A;i(ai)  and  kx{h2)  approach  the  single  crack  values.  As  shown  by  Boduroglu  and 
Erdogan  (1983)  and  as  may  be  observed  in  Figure  7,  for  a2  — >  hi  ki{ai)  and  ^1(62)  are 
highly  ill-defined  ( that  is,  their  slopes  are  unbounded  as  (ci  +  C2)/d  1). 

Some  sample  results  for  stress  intensity  factors  in  an  orthotropic  inhomogeneous 
medium  containing  two  unequal  collinear  cracks  under  plane  strain  conditions  are  given  in 
Figures  17-21  and  Tables  8  and  9.  Aside  from  the  geometric  parameters  C2/C1  and 
(ci  -f  C2)/d,  as  pointed  out  previously,  the  main  variables  in  this  case  are  the  material 
inhomogeneity  constant  a  and  the  shear  parameter  /c.  Both  the  figures  and  the  tables  show 
that  the  dependence  of  the  stress  intensity  factors  on  a  is  strong  and  on  k  is  rather  weak. 

The  results  for  the  crack/contact  problem  under  remote  loading  £22  =  £0  ei{x j a) 
are  shown  in  Tables  10,  11  and  Figures  22-25.  Figure  22  shows  the  crack  surface 
displacement  for  eo/^i  =  /5  =  0  (remote  bending)  obtained  by  ignoring  the  crack  surface 
interpenetration.  Thus,  on  the  compressive  side  of  the  external  load  the  solution  gives 
negative  displacement  which  is  physically  unacceptable.  Figures  23  and  24  show  the  crack 
surface  displacement  for  =  0  obtained  by  taking  into  account  the  smooth  contact  of  the 
surfaces.  Note  the  slight  dependence  of  the  crack  surface  displacement  on  the  sign  of  £i . 
Table  10  shows  the  stress  intensity  factors  and  the  size  of  the  contact  zone  again  for 
remote  loading  £22  =  {‘^\  /o)  and  for  various  values  of  the  material  inhomogeneity 

constant  a,  {E{xi)  =  EoQ-xpiaxi),  see  (4)).  The  contact  occurs  along  2:2  =  0, 
-a  <  xi  <  -b,  (-a  <  -b)  for  £1  >  0  and  along  X2  =  0,  b  <  Xi  <  a,  {b  <  a)  for 
£1  <  0.  The  table  also  shows  the  (theoretical)  stress  intensity  factors  calculated  by 
disregarding  the  crack  surface  penetration  (no  contact  case).  Note  that  the  presence  of 
contact  tends  to  magnify  the  stress  intensity  factor  at  the  crack  tip  on  the  tension  side  ( 
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ki{a)  for  £i  >  0  and  ki{-a)  for  £i  <  0),  Note  also  the  highly  significant  influence  of  the 
material  inhomogeneity  constant  on  the  stress  intensity  factors. 

Some  sample  results  for  the  combined  loading  "tension"  and  "bending", 
£22  =  £o  +  £^i(aJi/a),  {so/ei  =  P)  are  shown  in  Figure  25  and  Table  11.  Figure  25  shows 
the  crack  surface  displacement  for  various  values  of  p.  For  P>  Pc  and  the  contact  zone 
size  becomes  zero,  that  is  6  =  a,  P  =  Pc  corresponding  to  ki  (-a)  =  0.  Also,  6  <  a  for 
P  <  Pc.  Table  11  shows  the  crack  contact  zone  sizes  and  the  corresponding  stress 
intensity  factors.  If  the  material  is  homogeneous  and  isotropic  the  crack/contact  problem 
can  be  solved  in  closed  form  which  is  given  in  Appendix  A. 

Finally,  a  sample  result  giving  the  distribution  of  cr^y(x,0)  along  the  net  ligament 
6i  <  re  <  a2  is  shown  in  Fig.  26.  Note  that  the  solution  given  in  this  study  is  that  of  the 
perturbation  problem  which  must  be  added  to  the  results  obtained  from  the  uncracked 
medium  in  order  to  find  the  total  solution. 

5.  Conclusions 

In  the  case  of  two  collinear  cracks  with  crack  tips  at  (ai,fci)  and  (a2,fc2)  and 
ai  <  bi  <  a2  <  b2,  bi  -  ai  :$>  b2  -  02  (that  is,  (ai,  61)  being  the  dominant  crack ),  if  the 
medium  is  homogeneous  (isotropic  or  orthotropic)  then  under  mode  I  loading  generally 
ki(bi)  >  ki(ai)  >  ki{a2)  >  ki{b2),  meaning  that  fastest  crack  growth  would  take  place 
at  the  inner  tip  of  the  dominant  crack.  As  the  ligament  a2  —  &i  decreases  and  tends  to  zero 
ki  (a2)  as  well  as  ki  {bi)  becomes  unbounded.  However,  if  the  medium  is  inhomogeneous, 
the  values  of  the  stress  intensity  factors  may  be  influenced  quite  significantly  by  the 
relative  positions  of  the  dominant  and  the  secondary  cracks  with  respect  to  the  direction  of 
material  property  variation. 

After  introducing  the  "averaged"  orthotropy  parameters  E,  v,  k,  and  5o,  in  mode  I 
crack  problems  for  an  inhomogeneous  orthotropic  medium  E{p)  and  60  act  as  scaling 
constants  and  u  and  k  do  not  seem  to  have  a  significant  influence  on  the  stress  intensity 
factors.  In  crack  problems  involving  graded  materials  the  dominant  factor  remains  to  be 
the  material  inhomogeneity  constant. 

In  mode  I  and  mode  III  crack  problems  for  isotropic  graded  materials  described  by 
ix{x)  =  /roexp(ax)  the  kernels  of  the  integral  equations  can  be  evaluated  in  closed  form. 
Consequently,  for  these  materials  highly  accurate  benchmark  solutions  can  be  obtained.  In 
this  article  the  mode  I  stress  intensity  factors  in  such  materials  are  provided  for  tension  and 
bending  through  fixed  grips  as  well  as  third  degree  polynomial  tractions  acting  on  the 
crack  surfaces. 
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APPENDIX  A 


Crack  Closure  for  a  Homogeneous  Isotropic  Medium 


Let  the  medium  contain  a  single  crack  of  length  2a  along  the  x  axis  and  be  subjected 
to  crack  surface  tractions  axy{x,0)  —  0,  ayy{x,0)  =  p{x;),  —a<x<a.  By  defining 
dv{x,  -\-Q)/dx  =  g(x),  the  integral  equation  of  the  problem  may  be  expressed  as 


1 

r  g(s)ds 

vrj 

-a  S-X 

subject  to 

1 

r 

TT  j 

1  9{s)ds 

—a 

=  —p{x),  —a  <  X  <  a, 


A1 


A2 


where  Eq  =  E  for  plane  stress  and  Eq  ~  E / {1  -  for  plane  strain  conditions.  By 
substituting 


X  =  ax',  s  =  as',  g{s)  =  G(s'),  p{x)  =  P{x'). 
A1  and  A2  become 

7rJ_i  s'-x'  Eq  y_i 

We  now  assume  that  the  solution  of  (A4)  is  of  the  form 

AnTn{s') 


A3 


A4 


./2 


A5 


Considering  the  orthogonality  condition 

1,  m  =  n  =  0 


1  f'TMEM 

\/l  -  3^ 


V 

T^J-l 


dq  =  {  1/2,  m  =  n  >  1, 

0,  m^n, 


A6 


from  (A4b)  it  may  be  seen  that  Aq  =  0.  Also,  by  substituting  from  (A5)  into  (A4a)  and  by 
using  the  properties 


-f'— 

T^J-i  (q  - 


Tn{q) 


{q  -  r)y/l  -  q^ 


0,  n  —  0,  |r|  <  1 

dq={Un-i{r),  n>l,  |r|  <  1, 

Gn{r),  n  =  0,l,...,  Irl  >  1, 


A7 
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Tn{q)  =  cos(n0),  Un{q)  = 


cosd  =  q, 


A8 


Gnir)  = 


r|  (^“  (|r-|/r)Vr2  -  l) 

r  i/r2  -  1 


we  obtain 


A9 


X;  AnU„.i(x') 

m=l 


-C/O 


AlO 


As  an  example  consider  now  the  remote  "bending"  of  the  medium  for  which  we  have 


P{x')  =  -EoEix'  =  -EoeiUi{x')/2. 


All 


From  (AlO)  it  then  follows  that 

Ai=0,  A2  =  -£!■,  Am  =  0,  m  >  3. 
and  the  solution  becomes 

eiTiis')  «i(2s'"-l) 

VT^  vT^  ^ 

giving  the  stress  intensity  factors 

ki (a)  =  ^£^0^1  y/a,  h {-a)  =  -  ^EqSi i/a. 


A12 


A13 


A14 


Since  the  negative  stress  intensity  factor  (A14b)  implies  interpenetration  of  the  crack 
surfaces,  the  solution  given  by  (A13)  is  not  valid.  The  correct  solution  may  be  obtained  by 
taking  into  account  the  crack  surface  contact  near  the  end  x  —  —a.  Let  the  contact  region 
be  —a<x<—b  with  b  being  an  unknown  constant  to  be  determined  from  the 
smoothness  condition  ki(-b)  =  0.  The  integral  equation  (Al)  and  the  single-valuedness 
condition  (A2)  may  then  be  modified  as 


1 

~  i 

rg{s)ds 

=  J-p(x),  -b  <x  <a, 

A15 

TT  j 

-b  s-x 

il/Q 

1 

i 

[  g{s)  ds  = 

=  0, 

A16 

ttJ 

-6 

By  defining  the  following  normalizations 

s  =  cS  +  d,x  =  cX  +  d,c=  d  =  ^,9{s)  =  G{S),  p{x)  =  P{X),  A17 


(A  15)  becomes 
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Applying  again  remote  "bending"  we  have 


A18 


(X\ 

c  d 

[-)=-E^ex 

-x  +  - 

\a/ 

a  a 

=  P{X) 


A19 


The  integral  equation  (A18)  then  becomes 


f  «f)f  =  _2.. 

Li  s-x 

a 

A20 


a  a 


Again  assuming  the  solution  as  given  by  (A5)  it  can  be  shown  that 

OO 

Yi  =  -2£, 

m=l 

By  observing  that  Uq{X)  =  1  and  Ui  (X)  =  2X,  from  (A21)  it  follows  that 
d  c 

Ai=  -  2ei-,  A2  =  -si-,  Am  =  0,  m>2. 

a  a 

Thus,  the  solution  becomes 

G(S)  =  -T  P  „  Ns  +  c(l  -  2S2)] , 


A21 


A22 


A23 


giving 


ki  (a)  =  EqEi  a/c 


Za  —  b 
4a 

Zb  —  a 


A24 


k{-b)  = 


A25 


From  the  smooth  contact  condition  ki  (—6)  =  0,  the  unknown  constant  b  is  found  to  be 


,  1 
b  =  —a. 

3 

From  A24  and  A26  it  then  follows  that 

kiia)  ^  ^  or  =  0.5443. 


A26 


3  £^0^1  \/®  3 


A27 


Figure  A1  shows  the  crack  surface  displacements  obtained  from  (A13)  by  ignoring  the 
interpenetration  of  crack  surfaces  and  from  (A23)  by  taking  into  account  crack  closure 
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along  -a  <  X  <  -b.  From  Figure  A1  and  equation  (All)  it  is  seen  that  the  crack  closure 
near  the  end  x  —  -a  leads  to  magnification  of  stress  state  near  the  crack  tip  a;  =  a  and 
there  is  approximately  a  9%  increase  in  the  stress  intensity  factor  ki  (a)  over  the  nominal 
value  obtained  by  ignoring  crack  surface  contact. 


x/a 


Figure  A1  Crack  surface  displacements  in  a  homogeneous  isotropic 
medium  subjected  to  remote  bending  eyy  =  Si(x/a)  and  obtained  by 
ignoring  crack  closure  and  by  taking  it  into  account;  vq  =  2eia. 

After  solving  the  perturbation  problem  (A15)  the  stress  <Jyy(x,  0)  on  the  plane  of  the 
crack  may  be  obtained  by  using  (A7),  (A9),  (All)  and  (A18)  and  by  adding  the 
homogeneous  stress  ayy  =  EqEi  (x/a)  as  follows: 


0,  —b<x<a, 

^yy(^r  O)  ^  TP  f  C  cf 

’  ^[A,G,(X)  +  A2G2(X)]  +  Eos,  (  "  A  +  ^ 


a  a 


X  = 


X  —  d 


c  = 


a  -\-b 


d  = 


a  —  b 


c  2  ’  ”  2 

or  by  observing  that  b  =  a/3  it  may  be  shown  that 


00  <  a;  <  —b,  a  <  x  <  oo,  A28 
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EqSi 


1 


0,  |X|  <  1, 

2X  +  1 


+ 


|X1  >  1. 


A29 


Figure  A2  shows  the  stress  distribution  expressed  by  (A29).  Note  that,  as  expected, 
cryy{x,  0)  <  0  for  x  <  —b. 


-1  -b/a  1  x/a 


X 


Figure  A2  The  stress  ayy  in  the  plane  of  the  crack  in  a  homogeneous 
isotropic  medium  under  remote  "bending"  Syy  —  ei{x/ a). 


Let  us  now  assume  that  a  homogeneous  isotropic  medium  containing  a  crack  along 
-a  <  X  <  a,  y  =  0  is  subjected  to  remote  “tension”  Syy  =  eo  as  well  as  “bending” 
Syy  =  Si  (x/a).  The  input  function  in  the  integral  equation  (Al)  would  then  be 


p(x)  =  -Eo 


£q  +  Si 


(!)) 


=  —EoSi 


/5  + 


(!)]■ 


'-?■ 

Si 


A30 


giving  the  solution  as  follows: 

G{s')^  .  ^  \AiTi{s')  +  A2T2{s')\,  s' =  s/a,  A31 

V  1  —  5'^ 

Al  ==  —2eiP,  A2  "==  A32 
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h  (a) 
■£'o£i\/® 

fci(-a) 

Eo£i  y/a 


-{20  + 1), 


^(20-1), 


A3  3 

A34 


=  -a<^<a.  A35 

The  solution  given  by  (31)-(35)  is  valid  provided  ki{—a)  >  0.  From  (A34)  the  critical 
value  of  the  strain  ratio  0  =  Eq/ei  for  which  h  =  a  or  k\{— a)  =  0  is  found  to  be 


(£o/ei)c  —  0c  —  1/2. 


A3  6 


Thus,  crack  closure  would  take  place  if  /3  <  /3c  (or  if  ki{—a)  <  0  ).  In  this  case  the 
solution  is  obtained  from  (A15)  and  (A30)  by  assuming  that  6  <  a  and  by  following  the 
procedure  outlined  by  (A17)-(A27).  It  may  then  be  shown  that 


G{S)  =  -^2—[ArTi{S)  +  A2T2{S)],  S  =  {s-  d)/c, 

A2  =  -ei(^), 


6 

a 


5(1+4/?), 


c 

a 


3(1 + 0)1 


3(1+^) 


3/2 


d 

a 


|(l-2/3). 


h(-b) 

EoEiy/a 


=  0, 


v{x,  0) 
vq 


Vo  =  2eia,  x  =  cX  +  d, 


A3  7 

A3  8 

A3  9 

A40 

A41 

|X|  <  1.  A42 


Figure  A3  shows  the  crack  surface  displacement  v{x,0),  —b<x<a,  the  size  of  the 
contact  zone  b  and  the  normalized  stress  intensity  factors  for  various  values  of  the  strain 
ratio  P  =  Eo/ El. 
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b/a 

ki{-a) 

ki{a) 

Eo£i 

E(i£i^/a 
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1 
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1/3 

0.0 
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xja 

Figure  A3  The  crack  surface  displacement,  the  contact  zone  size  and  the 
normalized  stress  intensity  factors  in  a  cracked  isotropic  homogeneous 
medium  subjected  to  the  remote  loading  £yy  =  Sq  +  Si  (x/a),  vq  =  2£ia. 
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APPENDIX  B 


Evaluation  of  the  Kernel 


For  an  isotropic  inhomogeneous  medium  containing  a  series  of  collinear  cracks  along  the 
X  axis  the  kernel  of  the  integral  equation  (15)  may  be  expressed  as 


2 

R{cc)  =  r  / 


k  -\-ia 

+  4(A:^  +  ioik) 


Q^^dk  =  2  -  +  R2{uj)  ,  uj  =  s-x 


Defining  A  =  a/2  and  observing  that  k"^  +  iak  =  {k  +  iA)'^  +  A"^,  from  (Bl)  it  follows 
that 

R(uj)  ^  \  r  B2 

*  \J +  1)-^^  +  (A;  +  iA)^ 

Changing  the  variables  k  =  AV  and  assuming  A  >  0,  (B2)  becomes 

R{uS)  =  4  -=iK±2i_  B3 

*  Y  + 1)  +  +  if 

or,  by  letting  U  =  V-\-i,v-\-l  =  C‘^  and  B  =  Au  we  find 


B  /»oo+i 


U  +  i 
'C72  + 


e'^^du. 


By  substituting  U  =  -U'  and  observing  that 

r  (-g'  +  O  e-‘«g'(-di7')  =  - 
Joo-i  x/C^  +  U'^ 


(U'-i) 

'c^  + 


(B4)  may  be  written  as 


R{p))  =  Ae^ 


1  \  (U  ~  i)  ^-iBU 


ij  +  m 


1\  {U^i) 


iJ 


=  2Ae^3? 


{U  +  i) 


We  now  consider 

_  {U  + 1)  g  >  0  B7 

Ji  iy/C^  +  m 


Referring  to  Figure  B1  we  have 


/ 


(ii^e-v.=  r+/+p%/+r 

i  V  Ji  Jr  Ji{R+l)  J^  Ji(C-s) 


=  0, 


B8 


For  R-*oo  and  £-^0  it  can  be  shown  that  the  integrals  over  the  arcs  7  and  F  in  (B8)  are 
zero  and  from  (B8)  it  follows  that 


iyfC^TTP 


y+^^-By 


dy 


2/+  1 


dy. 


B9 


21 


=  I"*'  B-  >  0, 


Let  us  now  assume  that  B  <0,  B  =  -B'  and  define  the  function 

oo+i 

+  i72 

Referring  to  Figure  B2  we  have 

it+i  /*  /* — i{C-\-£) 

(R-1) 


B14 


/ 


(-a  + 1)  tB' 

c 


ix/c2T^ 


•^dz  =  +/+/  +/  +  /  =0’ 

Ji  Jr  J-i{R-i)  J'i  J-i{c-£) 


B15 


Again,  for  R-t’OO  and  e-+0  it  can  be  shovm  that  the  integrals  over  the  arcs  7  and  T  in 
B 15  are  zero  and  I2  becomes 


j^oo+i 


(U  +  i) 

i^CP  + 


dU 


=/■ 


ic 


if  Q-^^'^dZ 

J-iC  h 


iVC^ 


or 


3?- 


j'oo+i 


(U  +  i) 
i^/(PTTp 


-iB'U 


dU 


-I 


(-thA 


c 

_  f  °°  {—CY  +  1)  ^-b'cy 

~  Ji  Vy2-i 

=  -CKi{B'C)  +  Ko{B'C). 

Thus,  for  any  real  constant  B,  positive  or  negative,  we  find 
R{to)  =  2Ae®|c^ifi(|BlC)  +  Ko{\B\C)^  =  +  -^2(0;)), 


dY 


A.  =  Oi/2,  C  =  v^l  +  z^,  B  —  Aoj  —  ~ 


B17 

B18 

B19 
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APPENDIX  C 


The  Antiplane  Shear  Problem  for  a  Graded  Medium  Containing  a 

Crack 


It  can  be  shown  that  the  integral  equation  corresponding  to  an  infinite  inhomogeneous 
isotropic  medium  containing  a  crack  and  subjected  to  anti-plane  shear  loading  is  given  by 


1 

TT 


1 


S  —  X 


+  Ma{x,s) 


g{s)  ds 


2_ 

f^Q 


—a  <  X  <  a.  Cl 


,  ,  dw(x,0) 


B  =  ^{s-x)  C4 

where  Ko{z)  and  Ki{z)  are  the  modified  Bessel  fimctions  of  order  zero  and  one 
respectively.  It  is  again  assumed  that  the  shear  modulus  of  the  medium  is  approximated  by 
IJi(x)  =  fi,Qexp{ax).  Note  that  the  kernel  of  the  integral  equation  (C3)  is  nearly  identical  to 
that  given  by  (20)  for  the  in-plane  problem.  Table  Cl  shows  some  sample  results  obtained 
fi'om  (Cl)  and  their  comparison  with  that  obtained  by  delale  (1985)  numerically.  Since  the 
kernel  is  in  closed  form,  (Cl)  can  be  solved  within  any  desired  degree  of  accuracy. 


Table  Cl  The  normalized  mode  III  stress  intensity  factor  for  an 
inhomogeneous  isotropic  medium  subjected  to  uniform  crack  surface 
traction,  ayz{x,0)  =  —po,  —a<x<a. 


Delale  [1985] 

Present 

aa 

hia) 

ksi-a) 

hia) 

hi-a) 

Poy/a 

Po^/a 

TOa/o 

0.0 

1.0 

1.0 

1.0 

1.0 

0.1 

1.024 

0.973 

1.0228 

0.9731 

0.2 

1.045 

0.944 

1.0427 

0.9443 

0.3 

1.063 

0.914 

1.0605 

0,9149 

0.4 

1.080 

0.884 

1.0763 

0.8857 

0.5 

1.095 

0.855 

1.0906 

0.8570 

0.75 

1.127 

0.787 

1.1205 

0.7894 

1.00 

1.153 

0,726 

1.1438 

0.7291 
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Figure  3  Stress  intensity  factors  at  the  tips  ai  and  6iof  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy(x,0)  = -Eoeoexp(ax),  ko  =  Boeo^/^, 
Eo  =  Eo/(l  -  u^),  C2/C1  =  0.25,  =  0.3. 
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Figure  4  Stress  intensity  factors  at  the  tips  ai  and  61  of  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy{x,0)  = —Eo£oeKp{ax),  ko  =  Eoeo^yc[, 
Eo  =  Eo/{\  -  1/2),  (C2  -h  ci)/d  =  0.75,  z/  =  0.3. 
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Figure  7  Stress  intensity  factors  at  the  tip  ai  of  two  unequal  cracks  in  an  isotropic  FGM 
under  fixed-grip  condition,  ayy{x,0)  — -Eo£oexp{ax),  = 

Eo  —  Eo/{l  -  u'^),  aci  =  0.5,  u  =  0.3. 


aci 


Figure  8  Stress  intensity  factors  at  the  tips  a2  and  62  of  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy{x,0)  =  -£Jo£oexp(aa:),  ko  =  EoSoy/^, 
Eo  =  Eq/(1  -  C2/C1  =  0.1,  V  =  0.3. 


i^igure  9  Stress  intensity  factors  at  the  tips  ai  and  6iof  two  unequal  cracks  in  an  isotropic 
"GM  under  fixed-grip  condition,  cryy{x,0)  = -Eo£o^yip{ax),  ko  =  Eo£o^/^, 
Eo  =  Eo/(l  -  v\  C2/C1  =  0.25,  V  =  0.3. 
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Figure  10  Stress  intensity  factors  at  the  tips  oi  and  hi  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  ayy(x,0)  —  — £^o£oexp(ax),  ko  =  Eo£o^y^, 
Eo  =  Eo/(l  -  u^),  (C2  +  ci)/d2  =  0.75,  u  =  0.3. 
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Figure  1 1  Stress  intensity  factors  at  the  tips  ai  and  bi  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  ayy{x,0)  =  —  £io£'oexp(aar),  ko  — 

%  =  Eo/{l  -  v\  C2/C1  =  0.25,  V  =  0.3. 
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Figure  12  Stress  intensity  factors  at  the  tips  ai  and  61  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  ayy(x,0)  =  —Eo£o&xp(ax),  ko  =  EQ£o^/^, 
Eo  =  Eo /{I  -  v^),  ctci  =  0.5,  u  =  0.3. 
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Figure  13  Stress  intensity  factors  at  the  tips  a2  and  62  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  cryy(x,0)  =  —Eo€oexp(ax),  ko  =  Eo£o.y^, 
Eo  =  Eo/(l  -  z^2),  C2/C1  =  0.1,  u  =  0.3. 
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Figure  14  Stress  intensity  factors  at  the  tips  ai  and  61  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-load  condition,  cryy(x,  0)  =  -po,  ko  =  po^/^,  =  0-25, 
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Figure  15  Stress  intensity  factors  at  the  tips  ai  and  feiof  two  unequal  cracks  in  an  isotropic 
FGM under  fixed-load  condition,  ayy{x,0)  =  -po,  ko  —  Po^/^,  C2/C1  =  0.25,  p  =  0.3. 
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Figure  17  Stress  intensity  factors  at  the  tips  ai  and  fciof  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy(x,0)  = —eo-Eo/SQexp(axi), 
ko  =  {eoEo/6^)./^,  C2/C1  =  0.25,  k  =  0.5,  u  =  0.3. 
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Figure  18  Stress  intensity  factors  at  the  tips  ai  and  61  of  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy{x,0)  =  —£oEo/6Qexp{axi), 
ko  =  {eQEol6l).Jci,  (c2  -h  ci)/d  ^  0.75,  n  =  0.5,  v  =  0.3. 
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Figure  19  Stress  intensity  factors  at  the  tips  ai  and  biof  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy{x,0)  =  —Eoeo/6Qexp{axi), 
ko  =  (Eo£o/^o)\/^’  (^2  +  ci)/d  =  0.5,  C2/C1  =  0.25,  u  ^  0.3. 
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Figure  20  Stress  intensity  factors  at  the  tips  ai  and  61  of  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy{x,0)  = -Eo£o/6QQxp{axi), 
ko  =  (J5'o£o/5o)\/^j  =  0.5,  C2/C1  =  0.25,  u  —  0.3. 
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Figure  21  Stress  intensity  factors  at  the  tips  ai  and  61  of  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy{x,0)  = -Eoeo/6'^Qxp{axi), 
ko  =  Eo€ol6‘^y/ci,  aci  =  0.5,  u  =  0.3,  k  =  0.5. 
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Figure  22  Crack  surface  displacement  for  an  isotropic  inhomogeneous  medium  under 
fixed-grip  condition,  disregarding  interpenetration  of  crack  surfaces, 
ayy{x,Q)  =  -Eoei{xi/a)Qxp{axi),  vq  =  2eia,  Eq  ^  Eq/{1  -  v^),  v  =  0.3. 
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Figure  23  Crack  surface  displacement  for  an  isotropic  inhomogeneous  medium  under 
fixed-grip  condition  in  the  presence  of  smooth  contact,  (Jyy{x,  0)  =  -EqEi  {x/ a)exp(ax), 
vq  =  2eia,  Eo  =  Eo/{l  -  =  0.3,  £1  >  0. 
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Figure  24  Crack  surface  displacement  for  an  isotropic  inhomogeneous  medium  under 
fixed-grip  condition  in  the  presence  of  smooth  contact,  Oyyix,  0)  =  -EqSx  (a;/a)exp(aa:), 
vq  =  2sia,  Eq  =  Eo/{l  —  v'^),u  =  0.3, ei  <  0. 
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Figure25  Crack  surface  displacements  for  an  infinite  isotropic  inhomogeneous  medium 
under  fixed-grip  condition  :  Eyy  =  Eq  +Ei{x/a),  Eq/ei  =  p,  (3c  —  0.4485,  aa  —  0.5, 
=  0.3. 


x(ci 

Figure  26  The  stress  distribution  in  the  plane  of  the  crack  in  an  inhomogeneous  isotropic 
medium  with  two  unequal  colinear  cracks  under  the  remote  loading  c^ix,  0)  =  EoEoe°‘^, 
aci  =  0.25,  C2/C1  =  0.1,  (ci  +C2)/d  =  0.25,  ayy{x,0):  the  total  stress,  aPy{x,0)  :  the 
perturbation,  cr^{x,  0):  stress  at  infinity  (the  applied  stress),  =  a^y  +  a“. 
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Table  1  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  crack  surface  tractions  given  by  Eq.  21  {v  =  0.3  ). 


Plane  Strain 

ki{a) 

ki{-a) 

ki{a) 

A:i(-a) 

ki{a) 

h  {-a) 

ki{a) 

k  (-a) 

aa 

Po\/a 

PD\/a 

Pi  \/a 

Pi\f^ 

P2\/a 

P2 

P3 

P3  \/a 

0.00 

1.0 

1.0 

0.5 

-0.5 

0.5 

0.5 

0.375 

-0.375 

0.01 

1.0025 

0.5000 

-0.5000 

0.5006 

0.4994 

0.3750 

-0.3750 

0.10 

1.0238 

0.4998 

-0.4998 

0.5060 

0.4935 

0.3749 

-0.3749 

0.25 

1.0567 

0.9334 

0.4989 

-0.4986 

0.5142 

0.4833 

0.3745 

-0.3743 

0.50 

1.1062 

-0.4944 

0.5267 

0.4665 

mmmm 

-0.3722 

0.75 

1.1504 

-0.4878 

0.5379 

0.4505 

EEUlfl 

-0.3689 

isa 

mm 

-0.4793 

0.5483 

nitKa« 

-0.3646 

-0.4590 

0.5668 

-0.3542 

mi 

mm 

0.4660 

-0.4371 

0.5833 

-0.3429 

aa 

ki{a) 

fci(-a) 

ki{a) 

h  {-a) 

ki  (a) 

Po\/a 

Po^/a 

Pi\[a 

Pi \f^ 

P2\/a 

P2  \/a 

Psy/a 

0.00 

1.0 

1.0 

0.5 

-0.5 

0.5 

0.5 

0.375 

-0.375 

0.01 

1.0025 

0.9975 

0.5000 

-0.5000 

0.5006 

0.4994 

0.3750 

-0.3750 

0.10 

1.0235 

0.9737 

0.4998 

-0.4998 

0.5059 

0.4934 

-0.3749 

0.25 

1.0553 

0.9324 

0.4989 

-0.4985 

-0.3720 

0.50 

1.1019 

0.8640 

0.4962 

-0.4941 

0.5256 

EjEfciU 

-0.3720 

0.75 

1.1421 

0.8002 

0.4923 

-0.4870 

0.5359 

0,4494 

0.3712 

-0.3685 

1.00 

1.1774 

0.7427 

0.4879 

-0.4781 

msm 

mm 

mam 

1.50 

1.2369 

0.6473 

0,4780 

-0.4570 

mm 

mm 

0.3639 

mm 

2.00 

1.2862 

0.5746 

0.4680 

-0.4345 

0.5749 

0.3859 

0.3590 

-0.3415 

Table  2  The  normalized  stress  intensity  factors  for  an  inhomogeneous 
isotropic  medium  under  fixed  -  grip  condition  {  Eq  =  Eq  for  plane  stress, 
Eo  =  Eq / {l  —  for  plane  strain ,  u  =  0.3,  see  Eq.  30). 


m 

o 

11 

Plane  Strain 

Plane  Stress 

aa 

fci(a) 

fcl(-o) 

ki(a) 

ki{-a) 

fci(a) 

fcl(-o) 

fci(a) 

kii-a) 

SoEoy/a 

eoEoy/a 

eiEoy/a 

ei^oi/a 

eoEoy/a 

eoEoy/a 

eiEoy/a 

eiSc/a 

0.00 

1.0 

1.0 

0.5 

-0.5 

1.0 

1.0 

0.5 

-0.5 

0.01 

1.0075 

0.9925 

0.5050 

-0.4950 

1.0075 

0.9925 

0.5050 

-0.4950 

0.10 

1.0764 

0.9264 

0.5523 

-0.4523 

1.0761 

0.9261 

0.5523 

-0.4522 

0.25 

1.1986 

0.8230 

0.6402 

-0.3886 

1.1972 

0.8219 

0.6401 

-0.3885 

0.50 

1.4290 

0.6710 

0.8151 

-0.3010 

1.4245 

0.6683 

0.8146 

-0.3009 

0.75 

1.7029 

0.5434 

1.0326 

-0.2325 

1.6940 

0.5395 

1.0311 

-0.2326 

1.00 

2.0332 

0.4379 

1.3029 

-0.1794 

2.0189 

0.4331 

1.2997 

-0.1795 

1.50 

2.9294 

0.2811 

2.0570 

-0.1064 

2.9012 

0.2760 

2.0482 

-0.1067 

2.00 

4.2933 

0.1786 

3.2297 

-0.0631 

4.2454 

0.1741 

3.2107 

-0.0633 

Table  3  The  effect  of  Poison's  ratio  on  the  normalized  stress  intensity  factors 
for  an  inhomogeneous  isotropic  medium  (j&o  —  Eq  for  plane  stress, 
Eq  =  Eo/{l  —  for  plane  strain,  seeEqs.  28  and  30  ). 


Plane  Strain 

aa  =  0.5 

Fixed 

-Grip 

Fixed-Load 

fci(a) 

fcif-a) 

*i(a) 

ki{-a) 

ki{a) 

kii-a) 

ki{a) 

fci(-a) 

v 

eoEo^/^ 

e^EosJa 

€iEoy/a 

€lEo^/a 

Po\/a 

P\\/a 

Pi\/a 

0.0 

1,4129 

0.6615 

0.8131 

-0.3009 

1.0906 

0.8570 

0.4961 

-0.4932 

0.1 

jBsm 

-0.3009 

1.0950 

0.8597 

0.4961 

-0.4935 

-0.3009 

1.1001 

0.8629 

0.4962 

-0.4939 

0.3 

1.4290 

-0.3010 

1.1062 

0.8667 

0.4962 

-0.4944 

1.4368 

-0.3010 

1.1137 

0.8715 

0.4963 

-0.4950 

0.5 

1.4468 

0.6814 

0.8176 

-0.3011 

1.1234 

0.8775 

0.4964 

-0.4959 

Plane  Stress 

aa  =  0.5 

Fixed-Grip 

Fixed 

l-Load 

fcifa) 

ki{-a) 

All  (a) 

fci(-a) 

fei(a) 

/ci(-a) 

ki{a) 

ki{-a) 

V 

eQEo^/a 

€oEo^/a 

€lEo^/a 

Pi\/a 

Pi-v/a 

0.0 

1.4129 

0.6615 

0.8131 

-0.3009 

1.0906 

0.8570 

0.4961 

-0.4932 

0.1 

0.6639 

0.8136 

-0.3009 

1.0946 

0.8594 

0.4961 

-0.4935 

■ma 

0.4961 

-0.4938 

0.3 

1.4245 

■■niia 

0.4962 

-0.4941 

0.4 

Uligj 

TOaaiiiB 

0.4962 

-0.4943 

0.5 

D&EI 

-0.4946 

Plane  Strain 

Fixed-Grip 

Fixer 

-Load 

fci(a) 

fci(-a) 

fci(a) 

Alfa) 

Aif-a) 

Alfa) 

Aif-a) 

eoEo^^ 

elEo^/a 

e\Eo^/^ 

Po^/a 

Po\/a 

PiV“ 

Pi-^A 

0.0 

1.9819 

-0.1798 

1.1438 

0.7291 

0.4888 

-0.4752 

0.1 

1.9963 

0.4884 

-0.4763 

2.0132 

-0.1795 

1.1722 

0.7406 

0.4880 

-0.4777 

0.3 

2.0332 

-0.1794 

1.1902 

0.7480 

0.4876 

-0.4793 

0.4 

2.0576 

0.4459 

-0.1793 

1.2121 

0.7569 

0.4874 

-0.4813 

0.5 

2.0885 

0.4560 

ngn^H 

-0,1793 

1.2396 

0.7683 

0.4873 

-0.4840 

Plane  Stress 

aa  —  1.0 

Fixed-Grip 

Fixer 

-Load 

z/ 

ki{a) 

ki{-a) 

fci(a) 

ki{-a) 

Alfa) 

Aif-a) 

Poi/a 

Alfa) 

Pi\/a 

ki{-a) 

Pi\/a 

eoEoy/a 

eoEo^/^ 

^lEo^/a 

€iEoy/a 

pO'sj'a 

0.0 

0.4208 

-0.1798 

1.1438 

0,7291 

0.4888 

-0.4752 

0.1 

0.4251 

-0.1797 

1.1556 

0.7339 

0.4884 

-0.4762 

0.2 

mmmwm 

1.2971 

-0.1796 

_ 

1.1668 

0.7384 

0.4881 

1.2997 

1.1774 

0.7427 

0.4879 

0.4 

2.0301 

-0.1794 

1.1874 

0.7468 

0.4877 

-0.4790 

0.5 

2.0407 

0.4404 

-0.1793 

1.1970 

0.7507 

0.4875 

-0.4799 
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Table  4  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  fixed  -  load  condition  for  the  case  of  plane  strain  (z^  =  0.3  ). 


Present 

Delale  andErdogan  [1983] 

aa 

ki{a) 

po^/a 

ki{-a) 

Po^/a 

ki{a) 

piy/a 

ki{-a) 

Pi^/a 

ki{a) 

Po^/a 

ki{-a) 

Po\A 

ki  (a) 

pi^a 

k  {-a) 

Pl^/a 

0.00 

1.0 

1.0 

0.5 

-0.5 

1.0 

1.0 

0.5 

-0.5 

0.01 

1.0025 

0.9975 

-0.5000 

1.003 

0.10 

-0.4998 

1.026 

E^HI 

SQQIII 

0.25 

0.4989 

-0.4986 

1.061 

0.931 

0.498 

0.8667 

-0.4944 

0.863 

0.494 

■gna 

ussm 

0.8044 

0.801 

0.489 

0.7480 

bee^i 

00911 

0.483 

-0.481 

Table  5  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  fixed  -  load  condition  for  the  case  of  plane  stress  (i/  =  0.3  ). 


Present 

Delale  and  Erdogan  [1983] 

ki{a) 

Po\/a 

A:i(-a) 

Po^/a 

ki{a) 

Pi^/a 

ki{-a) 

Pi^/a 

ki{a) 

Po^/a 

ki{-a) 

po^/a 

h  (a) 

Pi'/a 

k  (-a) 

Pi\/a 

1.0 

1.0 

0.5 

-0.5 

1.0 

-0.5 

BliW 

mi 

mm 

0.500 

mmm 

QQ 

0.9737 

0.4998 

mm 

mm 

1.0553 

0.9324 

0.4989 

1.060 

1.1019 

0.8640 

■igfeRfcM 

beeeo 

1.113 

ibesi 

0.75 

0.8002 

BE^II 

0.797 

-0.489 

1.00 

1.1774 

0.7427 

0.4879 

-0,4781 

1.209 

0.740 

0.483 

-0.480 

Table  6  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  fixed  -  grip  condition  (  Eq  =  Eq  for  plane  stress, 
=  Eq /{I  —  for  plane  strain  ). 


Plane  Strain 

Plane  Stress 

o 

II 

Present 

Delale  and  Erdogan  [1983] 

Present 

Delale  and  Erdogan  [1983] 

aa 

ki{a) 

fci(-a) 

ki{a) 

ki{  a) 

ki(a) 

hi-a) 

ki{a) 

ki(-a) 

eoEo^/a 

e^E^y/a 

eoEoy/a 

e^Eoy/a 

eoEoy/a 

eoEoy/a 

eoEoy/a 

eoEo^ 

0.00 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.01 

1.0075 

0.9925 

1.008 

0.992 

1.0075 

0.9925 

1.008 

0.10 

1.078 

0.925 

■lilrMM 

0.9261 

1.078 

0.25 

BEn 

1.203 

0.821 

0.8219 

1.202 

mm 

0.50 

1.4290 

0.6710 

1.439 

0.667 

1.4245 

0.6683 

1.435 

0.665 

0.75 

1.7029 

0.5434 

1.721 

0.539 

1.6940 

0.5395 

1.713 

0.535 

1.00 

2.0332 

0.4379 

2.063 

0.433 

2.0189 

0.4331 

2.048 

0.429 
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Table  7  Stress  intensity  factors  at  the  tips  x  =  ax  and  x  —  h\  of  two  unequal 
cracks  in  an  isotropic  FGM  for  the  case  of  plane  strain,  jy  =  0.3, 
Eo  =  Exx/{1  -  v'^). 


Fixed-Grip 

(yyy{x,  0)  =  -.EoS^oexp(Q!a;),  Hq  =  EoEq^i 

aci  =  0.5 

C2/C1  =  0.25 

C2/C1  =  0.75 

C2/C1  =  1.0 

(cj  +  C2)/d 

ko 

fci(fci) 

ko 

fci(ai) 

ko 

kiih) 

ko 

ki{ai} 

ko 

kiih) 

ko 

0.0 

0.6710 

1.4290 

0.6710 

1.4290 

0.6710 

1.4290 

0.1 

0.6713 

1.4299 

0.6728 

1.4334 

0.6735 

1.4350 

0.2 

0.6720 

1.4319 

0.6766 

1.4435 

0.6791 

1.4495 

0.4 

0.6740 

1.4405 

0.6873 

1.4805 

0.6947 

1.5000 

0.6 

0.6773 

0.7030 

1.5601 

0.7171 

1.6032 

0.8 

0.6835 

■fsEgiH 

0.7284 

1.7721 

0.7526 

1.8624 

0.975 

0.7034 

2.3181 

0.7883 

3.1421 

0.8314 

3.4198 

1.0 

0.7656 

00 

0.9395 

00 

1.0210 

00 

Fixed-Load 

^1/2/ (^jO)  Po,  ^0  Pos,/^ 

ctci  =  0.5 

(ci  +C2)/d 

0.0 

0.8667 

1.1062 

0.8667 

1.1062 

0.8667 

1.1062 

0.1 

0.8667 

0.8667 

0.8667 

1.1062 

0.2 

0.8668 

0.8668 

■miMW 

0.8668 

1.1064 

0.4 

0.8674 

0.8687 

mmvem 

0.8688 

1.1125 

0.6 

0.8692 

bieb 

0.8750 

miB 

0.8765 

1.1436 

0.8 

0.8735 

0.8893 

1.2445 

0.8945 

1.2584 

0.975 

0.8881 

BsiS&Si 

0.9280 

2,0652 

0.9422 

2.1280 

1.0 

0.9334 

00 

1.0258 

00 

1.0503 

00 

Table  8  The  normalized  stress  intensity  factors  for  two  unequal  cracks  in  an 
infinite  inhomogeneous  orthotropic  medium  under  uniform  crack  surface 
pressure,<7y2/(a;, 0)  =  -pQ,E{xi)  =  Eoex^{axi),v  =  0.3,  (ci  +  C2)/d  =  0.5. 


K  =  -0.25,  C2/C1  =  0.25 

K  =  0.0,  C2/C1  =  0.25 

a  Cl 

ki{ai) 

Po^/^ 

kiih) 

Po^/^ 

kiiao) 

Po^/^ 

kiih) 

POy^ 

aci 

ki  (ai) 
Po^/^ 

h{bi) 

po-v/cT 

hiao) 

Poa/cT 

klih) 

Po^/^ 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

0.01 

1.00171 

1.01253 

0.55240 

0.54232 

0.01 

1.00171 

1.01252 

0.55240 

0.54231 

0.10 

0.97825 

1.03400 

0.55437 

0.54991 

0.10 

0.97786 

1.03357 

0.55421 

0.54976 

0.25 

0-93900 

1.06872 

0.55710 

0.56227 

0.25 

0.93743 

1.06671 

0.55646 

0.56166 

0.50 

0.87545 

1.12420 

0.56016 

0.58189 

0.50 

0.87163 

1.11809 

0.55856 

0.58038 

0.75 

0.81612 

1.17694 

0.56131 

0.60011 

0.75 

0.81038 

1.16556 

0.55884 

0.59774 

1.00 

0.76203 

1.22710 

0.56061 

0.61679 

1.00 

0.75489 

1.20973 

0.55749 

0.61377 

1.25 

0.71342 

1.27482 

0.55823 

0.63194 

1.25 

0.70538 

1.25104 

0.55471 

0.62849 

1.50 

0.67015 

1.32028 

0.55437 

0.64562 

1.50 

0.66160 

1.28987 

0.55069 

0.64197 

1.75 

0.63181 

1.36367 

0.54923 

0.65794 

1.75 

0.62305 

1.32655 

0.54562 

0.65430 

2.00 

0.59792 

1.40519 

0.54304 

0.66901 

2.00 

0.58915 

1.36137 

0.53969 

0.66557 
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Table  8  (Continued ) 


A€  =  0.5,  C2/C1  =  0.25 

/£  =  1.0,  C2/C1  =  0.25 

aci 

kM 

aci 

fci(ai) 

kiih) 

ki{a2) 

kM 

Po^ 

Po^/^ 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

UjO^yl 

1.01010 

0.55216 

0.54146 

0.01 

1.00170 

1.01251 

0.55239 

0.54231 

lim 

1.01251 

0.55239 

0.54231 

0.10 

0.97746 

1.03312 

0.55404 

0.54960 

0.10 

0.97725 

1.03289 

0.55395 

0  .54952 

0.25 

0.93578 

1.06459 

0.55575 

0.56099 

0.25 

0.93491 

1.06348 

0.55536 

0.56061 

ii£ii 

0.86754 

1.11153 

0.55669 

0.57857 

0.50 

0.86538 

1.10805 

0.55563 

0.57753 

0.75 

0.80420 

1.15325 

0.55579 

0.59474 

0.75 

0.80091 

1.14668 

0.55399 

0.59294 

1.00 

0.74718 

1.19085 

0.55341 

0.60968 

Hjjjl 

1.18072 

0.55090 

0.60712 

1.25 

0.69669 

1.22509 

0.54981 

0.62349 

Og 

1.21113 

0.54669 

0.62021 

1.50 

0.65236 

1.25657 

0.54521 

0.63626 

Ijgi 

0.64742 

1.23862 

0.54157 

0.63234 

1.75 

0.61357 

1.28579 

0.53981 

0.64809 

ggi 

0.60850 

1.26377 

0.53577 

0.64362 

2.00 

0.57964 

1.31314 

0.53375 

0.65906 

2.00 

0.57457 

1.28702 

0,52942 

0.65412 

K  =  5.0,  C2/C1  =  0.25 

K,  =  —0.25,  C2/C1  —  1.0 

aci 

A:i(ai) 

Po^/^ 

ki{bi) 

Po^/^ 

fci(a2) 

aci 

fci(ai) 

Po^/^ 

ki{a2) 

POy/^ 

^1(^2) 

POy/^ 

0.0 

1.00429 

1.01010 

0.55216 

0.0 

1.02796 

1.04795 

1,02796 

0.01 

1.00169 

1.01250 

0.55239 

HQI 

1.02477 

1.04978 

1.0311 

0.10 

0.97681 

1.03240 

0.55376 

IfUil 

.99602 

1.06556 

li^a 

1.0591 

0.25 

0.93307 

1.06111 

0.55451 

0.55979 

IQ 

.94961 

1.09128 

0.99839 

0.50 

0.86075 

1.10060 

0.55317 

0.57510 

.87927 

1.13545 

0.94573 

lUU 

0.79385 

1.13250 

0.54967 

0.58855 

liQ 

.81729 

1,18161 

loa 

0.73421 

1.15878 

0.54467 

0.60059 

UjQI 

.76231 

1.22847 

laa 

mwjM 

0.68203 

1.18078 

0.53860 

0.61150 

1.25 

IQQ 

loa 

mil 

0.63676 

1.19948 

0.53179 

0.62148 

1.50 

QjQI 

■iWaiTSM 

ima 

IQ 

0.59757 

1.21562 

0.52447 

0.63068 

1.75 

.63179 

1.36299 

0.71025 

usm 

0.56361 

1.22978 

0.51680 

0.63921 

2.00 

.59791 

1.40459 

0.67350 

K  =  0.0,  C2/C1  =  1.0 

K  =  0.5,  C2/C1  —  1.0 

aci 

ki{ai) 

fci(6i) 

fci(a2) 

h{b2) 

aci 

fc](ai) 

kiibi) 

^1(02) 

kiib2) 

POy/^ 

Poy^ 

Po^/^ 

POy^ 

POy/^ 

Po^/^ 

Poy^ 

Poy^ 

0.0 

1.02796 

1.04795 

1.04795 

1.02796 

0.0 

1.02796 

1.04795 

1.04795 

1.02796 

0.01 

1.02476 

1.04977 

1.0461 

1.0311 

0.01 

1.02475 

1.04975 

1.0461 

1.0311 

0.10 

0.99552 

1.06498 

1.0282 

1.0585 

0.10 

0.99499 

1.06435 

1.0276 

1.0579 

0.25 

0.94786 

1.08894 

0.99606 

1.0190? 

0.25 

0.94598 

1.08641 

0.99350 

1.0990 

0.50 

0.87537 

1.12905 

0.94040 

1.1691 

0.50 

0.87115 

1.12210 

0.93434 

1.1608 

0.75 

1.2306 

0.75 

0.80530 

1.15758 

0.87674 

1.2159 

1.00 

1.2865 

1.00 

0.74746 

1.19212 

0.82307 

1.2651 

1.25 

0.70543 

1.25099 

0.78567 

1.25 

0.69673 

1.22504 

0.77421 

1.3092 

1.50 

0.66159 

1.28932 

0.74185 

1.50 

0.65235 

1.25606 

0.73029 

1.3489 

1.75 

0.62304 

1.32590 

0.70217 

1.4255 

1.75 

0.61356 

1.28518 

0.69109 

1.3848 

2.00 

0.58914 

1.36080 

0.66641 

1.4640 

2.00 

0.57964 

1.31259 

0.65616 

1.4175 
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Table  8  (Continued ) 


K  =  1.0,  C2/C1  =  1.0 

K  =  5.0,  C2/C1  =  1.0 

fci(ai) 

ki{a2) 

k,{b2) 

POy/^ 

a  Cl 

WtfiM 

1.02796 

1.04795 

1.04795 

1.02796 

0.0 

tmm 

0.01 

1.02474 

1.04975 

1.0461 

1.0311 

0.01 

1.02473 

1.04973 

1.0461 

1.0311 

0.10 

0.99471 

1.06402 

1.0272 

1.0575 

0.10 

0.99412 

1.06333 

1.0265 

1.0568 

0.25 

0.94499 

1.08506 

0.99213 

1.0974 

0.25 

0.94284 

1.08218 

0.98913 

1.0940 

0.50 

0.86890 

1.11838 

0.93097 

1.1562 

0.50 

0.86404 

1.11032 

0.92338 

1.1457 

0.75 

0.80198 

1.15086 

0.87168 

1.2075 

UggI 

0.79481 

1.13628 

0.85991 

1.1881 

1.00 

0.74334 

1.18194 

0.81682 

1.2527 

0.73445 

1.15979 

0,80183 

1.2234 

1.25 

0.69209 

1.21107 

0.76727 

1,2927 

1.25 

0.68207 

1.18065 

0.75005 

1.2530 

1.50 

0  .64741 

1.23812 

0.72308 

1.3285 

1.50 

0.63675 

1.19898 

0.70449 

1.2784 

1.75 

0.60849 

1.26318 

0,68390 

1.3606 

1.75 

0.59756 

1.21507 

0.66462 

1.3002 

2.00 

0.57456 

1.28649 

0.64923 

1.3897 

2.00 

0.56360 

1.22929 

0.62976 

1.3194 

Table  9  The  normalized  stress  intensity  factors  for  two  unequal  cracks  in  an 
infinite  inhomogeneous  orthotropic  medium  under  uniform  crack  surface 
pressure,  (Tyy(x,  0)  =  -po,  E{xi)  =  EoeKp{axi),  v  =  0.3,  k  =  0.5. 


C2/C1  —  0.75,  (ci  4-  C2)/d  —  0.75 

kiih) 

POy/^ 

fci(a2) 

POy/^ 

h(h) 

Poy^ 

0.0 

1.14076 

1.04573 

0.94551 

0.01 

1.14252 

1.0452 

0.94823 

0.10 

1.15652 

1.0389 

0.97155 

0.25 

0.96426 

■raaaa 

1.0244 

1.0077 

0.50 

Em 

1.0632 

0.75 

0.81314 

1.22513 

1.1138 

■BBl 

0.75192 

1.24726 

0.92738 

1.1604 

0.69914 

1.26906 

0.89413 

■Kbl 

0.65361 

1.29073 

0.86230 

eeebi 

1.75 

0.61421 

1.31230 

0.83227 

1.2802 

2.00 

0.57996 

1.33373 

0.80422 

1.3146 

C2/C1  =  0.25,  (ci  +  C2)/d  =  0.75 

kiih) 

P0\^ 

0.0 

EEQEa 

1.04850 

0.66867 

0.60894 

0.01 

1.05090 

0.66935 

0.61018 

0.10 

1.07123 

0.67479 

0.62088 

0.94123 

1.10193 

EKI 

nKia 

0.87140 

E^a 

0.75 

0.80678 

0.69652 

0.68704 

EEB 

0.74884 

1.22084 

0.69956 

0.70907 

EE&I 

0.69772 

1.25217 

0.70055 

0.72964 

0.65298 

1.28079 

0.69990 

0.74892 

1.75 

0.61394 

1.30729 

0.69792 

0.76707 

2.00 

0.57986 

1.33211 

0.69487 

0.78420 

C2/C1 

=  1.0,  ( 

Cl  +  C2)/d  =  0.75 

aci 

ki (ai) 

fci(fei) 

ki  (02) 

ki{b2) 

Po^/^ 

Po^/^ 

Po^/^ 

0.0 

1.06844 

1.17318 

1.17318 

1.06844 

0.01 

1.06475 

1.17451 

1.1717 

1.0720 

0.10 

1.03045 

1.18417 

1.1565 

1.1028 

0.25 

0.97370 

1.19589 

1.1262 

1.1500 

0.50 

0.88781 

1.21176 

1.0707 

1.2213 

0.75 

0.81450 

1.22790 

1.0150 

1.2853 

1.00 

0.75226 

1.24581 

0.96240 

1.3434 

1.25 

0.69915 

1,26541 

0.91427 

1.3963 

1.50 

0.65354 

1.28619 

0.87088 

1.4448 

1.75  ! 

0.61414 

1.30762 

0.83206 

1.4894 

2.00 

0.57991 

1.32931 

0.79742 

1.5307 
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Table  9  (  Continued  ) 


C2/C1  =  0.25,  aci  =0.5 

C2/C1  =  0.75  ,  aci  =  0.5 

Cl  +  C2 

d 

ki(ai) 

ki(bi) 

POy^ 

ki{a2) 

kl{b2) 

Poi/cT 

Cl  +  C2 

d 

hih) 

Po^/^ 

ki{a2) 

po^y^ 

0.0 

0.8661 

1.1053 

0.4331 

0.5527 

0.0 

1.1053 

0.7501 

0.1 

0.86614 

1.10527 

0.48805 

0.51930 

0.1 

0.86614 

1.10526 

0.78301 

0.94138 

0.2 

0.86619 

1.10539 

0.49731 

0.52839 

0.2 

0.86622 

1.10545 

0.79236 

0.95090 

0.86640 

1.10607 

0.51092 

0.54096 

0.3 

0.86674 

1.10694 

0.80597 

0.96313 

0.4 

0.86683 

1.10783 

0.53004 

0.55743 

0.4 

0.86809 

1.11143 

0.82517 

0.97818 

0.5 

0.86754 

1.11153 

0.55669 

0.57857 

mm 

0.87055 

1.12107 

0.85224 

0.99643 

0.6 

0.86863 

1.11885 

0.59451 

0.60569 

0.87439 

1.13943 

0.89137 

1.0187 

0.7 

0.87026 

1.13374 

0.65093 

0.64115 

0.7 

0.88008 

1.17380 

0.95114 

1.0464 

0.8 

0.87287 

1.16752 

0.74431 

0.68972 

0.8 

0.88862 

1.24313 

1.0531 

0.9 

0.87776 

1.26963 

0.94244 

0.76440 

0.9 

0.90295 

1.42220 

1.2792 

1.1376 

0.95 

0.88273 

1.44262 

1.1969 

0.82657 

0.95 

0.91588 

1.1836 

1.0 

0.9326 

00 

00 

1.2606 

1.0 

1.0244 

00 

00 

1.5458 

1  C2/C1  =  1.0,  ac\  =  0.5 

ki  (ai) 
Voy/Ei 

fci(6i) 

PDa/^ 

ki  (a2) 

ki{b2) 

PO-y/Q 

0.8661 

1.1053 

0.8661 

1.1053 

0.86614 

1.10526 

0.86986 

1.1096 

0.86620 

1.10541 

0.87854 

1.1187 

0.86671 

1.10682 

0.89116 

1.1302 

0.86823 

1.11153 

0.90900 

1.1441 

0.87115 

1.12210 

0.93434 

1.1608 

0.87588 

1.14252 

0.97136 

1.1810 

0.88301 

1.18072 

1.0287 

1.2061 

0.89377 

1.25695 

1.1282 

1.2392 

0.91169 

1.45009 

1.3533 

1.2894 

0.92762 

1.73266 

1.6602 

1.3323 

1.0561 

00 

00 

1.6792 

Table  10  Stress  intensity  factors  and  the  contact  zone  sizes  for  an  isotropic 
inhomogeneous  medium  under  fixed  grip  remote  bending,  u  =  0.3, 

ko  =  %el^/a,  ayy{x,Qi)  =  -Eoffi (a:/a)exp(Q;x), Eo  =  Eo/{l  -  p^). 


El  >  0 

El  <  0 

No  contact 

aa 

hja 

ki{a) 

hja 

fci(-a) 

fci(a) 

fci(-a) 

EoEiv/o 

E{)e\yJ~a 

0.0 

0.3333 

0.5443 

0.3333 

0.5443 

0.5 

-0.5 

0.1 

0.3405 

0.5964 

0.3258 

0.4965 

0.5523 

-0.4523 

0.25 

0.3508 

0.6832 

0.3140 

0.4320 

0.6402 

-0.3886 

0.5 

0.3669 

0.8554 

0.2947 

0.3419 

0.8151 

-0.3010 

1.0 

0.3968 

1.3359 

0.2585 

0.2130 

1.3029 

-0.1794 

1.5 

0.4234 

2.0827 

0.2269 

0.1320 

2.0570 

-0.1064 

2.0 

0.4476 

3.2488 

0,2000 

0.0815 

3.2297 

-0.0631 
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Table  1 1  Stress  intensity  factors  and  contact  zone  sizes  for  an  isotropic  inhomogeneous 
medium  under  remote  strain  Syy  =  £o  +  =  £i{P  +  a:/a), (Fig.  25),  a  =  0.5, 

u  =  0.3,  eo/ei  =  (3,  £i>  0. 


^=7 

hja 

k^{a) 

ki{-a) 

Eo£\y/a 

Eo£iy/a 

0.0 

0.3669 

0.8554 

0.0 

0.25 

0.7184 

1.1798 

0.0 

0.4485 

1 

1.4561 

0.0 

0.45 

1 

1.4582 

0.0010 

0.50 

1 

1.5297 

0.0345 

0.625 

1 

1.7083 

0.1184 

0.75 

1 

1.8869 

0.2023 

1.00 

1 

2.2442 

0.3700 
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